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摘 要 : 从 流域 尺度 揭示 天 山北 坡 六 川 变化 状况 ,对 下 游 绿 洲 地 区 水 资源 的 合理 开发 利用 具有 重要 
意义 。 基 于 多 源 遥 感 影像 提取 天 山北 坡 诸 河流 域 近 期 冰川 这 界 ,结合 前 期 发 布 的 冰川 编目 及 AS- 
TER DEM 数据 ,对 该 区 域 冰川 面积 和 物质 平衡 变化 特征 进行 了 分 析 。 结 果 表 明 :(1) 2015 年 前 后 


天 山北 坡 共计 分 布 冰川 10061 条 ,总 
xk Jl i8 
1999 一 201$ 年 天 山北 坡 国外 段 冰川 
山北 坡 冰川 表 
0.04 m w.e. a (3I) 天 山北 坡 东 、 


Bol 


积 约 4855.85+245.86 km, 1960s—2015 年 天 山北 坡 国 内 段 
45 2E 4i SE HE 7g 0.5296 - 710.0690 a^^, EL XE 4E R E dr x XR a 3 (0.9696 + a'+0.88% a!) ; 
P E Ag ik R 21 7] 0.56%- a^ 30.3196- a^, (2) 2000—2020 年 天 
高 程 变化 速率 约 为 -0.57+0.01 m'.a ,冰川 物质 呈 持 续 亏 损 态 ,物质 平衡 为 -0.39+ 
只 和 物质 平衡 变化 均 存在 一 定 的 空间 差异 性 , 近 十 儿 


年 间 , 东 段 各 子 流 域 冰川 面积 萎缩 速率 和 物质 亏损 速率 均 相对 较 大 。 


关 键 词 : 天 山北 坡 ; 冰川 变化 ; 物质 平衡 ;遥感 
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冰川 作为 冰冻 圈 的 重要 组 成 部 分 ,对 气候 变化 
的 响应 极为 敏感 1。IPCC (Intergovernmental Panel 
on Climate Change) 第 六 次 评估 报告 (AR6) 指 出 ,最 
近 10 a 间 (2011 一 2020 年 ) 地 球 表层 平均 气温 比 
1850—1900 年 高 出 1.09 % (0.95~1.20 €), HIE 
50a(1970s 以 来 ) 是 近 2000 a 来 气温 升幅 最 快 的 50 a?! o 
在 气候 变 暖 的 背景 下 ,全 球 山 地 冰川 普遍 呈 退 缩 趋 
势 , 造 成 海平 面 上 升 .极端 气候 事件 增多 .灾害 频 发 
等 一 系列 气候 环境 问题 ,引起 了 国际 社会 的 广泛 关 
注 汪 5。 此 外 ,全 球 约 7$% 的 淡水 资源 都 来 自 于 冰 
川 ,冰川 变化 会 对 区 域 水 资源 平衡 和 经 济 社会 可 持 
续 发 展 带 来 重要 影响 。 因 此 ,冰川 变化 研究 不 仅 
具有 重要 的 科学 意义 ,同时 也 对 区 域 经 济 社会 发 展 
具有 重大 的 现实 意义 。 

天 山 山脉 地 处 中 亚 干 旱 区 ,是 世界 上 山地 冰川 
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分 布 最 广 的 山系 之 一 。 冰 川 融 水 是 该 区 域 重 要 的 
水 资源 补给 ,气候 变 暧 导致 的 冰川 加 速 消融 和 退缩 
对 区 域 水 资源 量 及 年 内 分 配 产生 了 重要 影响 ””。 
因此 ,系统 开展 该 区 域 冰川 变化 研究 ,对 下 游 绿 洲 
地 区 未 来 水 资源 的 合理 规划 、 利 用 与 管理 具有 重要 
指导 意义 ,还 能 为 促进 区 域 可 持续 发 展 提供 重要 科 
学 依据 。 目 前 ,国内 外 专家 学 者 已 经 对 天 山地 区 的 
冰川 变化 进行 了 深入 的 研究 。 邢 武 成 基于 Land- 
satTMETM+ 遥 感 影像 提取 2 期 中 国 天 山 冰川 边界 ， 
分 析 冰 川 变 化 特征 ,研究 发 现 1989 一 2017 年 中 国 天 
山 冰 川 总 体 呈 萎缩 状态 ,冰川 面积 减少 约 879.73 km’, 
蓉 缩 速率 达到 0.52%.a, 且 东 天 山 冰川 面积 萎缩 速 
率 最 快 。Chen 等 ”对 天 山 40 个 流域 的 冰川 变化 数 
据 进 行 统计 分 析 ,发现 1960s 一 2010s 天 山 约 97.52% 
的 冰川 星 退 缩 状 态 ,其 中 天 山北 部 和 东部 退缩 趋势 
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尤为 明显 。Farinott 等 "基于 重力 卫星 数据 .激光 
测 高 数据 以 及 天 山 不 同 区 域 的 冰川 面积 数据 ,对 整 
个 天 山 的 冰川 物质 亏损 状况 进行 评估 ,研究 发 现 
1961 一 2012 年 天 山 冰 川 总 面积 减少 约 18%+6%, 物 
质 亏 损 约 27%+15%。 此 外 ,还 有 众多 学 者 对 天 山 不 
同 区 域 或 子 流域 如 东 天 山中 、 伊 犁 河流 域 5 Ter- 
skey- Alatoo 区 域 " 等 的 冰川 变化 进行 研究 , 均 发 现 
研究 区 内 冰川 普遍 呈 退 缩 趋势 。 

已 有 研究 主要 侧重 于 天 山 国内 段 不 同 子 流域 
的 冰川 变化 特征 。 本 文 依据 河流 走向 ,以 流域 尺度 
划分 天 山北 坡 范 围 , 基 于 多 源 遥 感 资料 获取 天 山北 
坡 诸 河流 域 近期 的 冰川 边界 ,结合 前 期 发 布 的 冰川 
编目 数据 ,揭示 该 区 域 近 60 a 来 冰川 面积 时 空 变化 
特征 ,并 基于 ASTER DEM 数据 分 析 近 20 a 来 该 区 
域 冰 川 物质 平衡 时 空 变化 特征 ,以 期 对 整个 天 山北 
坡 冰 川 变 化 状况 有 更 全 面 的 认识 ,从 而 为 下 游 绿 洲 
地 区 水 资源 的 合理 开发 利用 提供 重要 的 理论 文 撑 。 


1 研究 区 概况 


天 山 山脉 位 于 亚洲 内 陆 腹 地 ,是 亚洲 中 部 最 大 

i 也 是 世界 上 山地 冰川 分 布 最 多 的 山系 
。 天 山 山 脉 西 起 图 兰 平原 ,向 东 穿 越 吉 尔 吉 其 

ind 国 新 疆 境内 ,全 长 约 2100 km, 


zf Kara- Batkaky Me 
图 例 T^ j 
© 主要 城市 ”~ 一 河流 x 监测 冰川 
es 冰川 于 湖泊 ”一 流域 边界 


一 天 山北 坡 界线 0 — 100km 


注 :图 中 数字 代表 流域 编码 。1 : 伊 吾 河 流域 ;2 : 幻 彩 湖 流域 ;3: 巴 里 坤 湖 流域 ;4: 博 格 达 山 北 坡 流域 ;5: 


南北 最 大 宽度 约 300 km ^^ ,其 中 ,天 山 国 内 段 全 长 
约 1700 km, 山体 总 面积 约 5.7x10 km? LLL 
区 地 形 高 差 悬 殊 L^ 温 变 化 大 且 具 有 胃 显 的 空 x [8] 25 
5r , 北 坡 平 均 气温 低 于 南 坡 ;降水 主要 受 西风 气流 
和 北冰洋 气 团 的 影响 ,空间 分 布 不 均 , 北 坡 降水 量 
多 于 南 坡 ” 。 本 文 研究 区 为 天 山北 坡 (图 1) ,主要 包 
括 哈 尔 里 克 山 、 巴 里 坤 山 、 博 格 达 山 、 天 格 尔 山 、 依 
连 哈 比 尔 藉 山 、 那 拉 提 山 、 科 克 铁 克 山 、 哈 尔 克 他 马 
山 \ 捷 尔 斯 格 伊 阿拉 套 和 塔 拉 斯 阿拉 套 的 北 坡 , 以 
及 博 罗 克 努 山 、 科 古人 鞭 山 、 别 珍 套 山 、 外 伊犁 阿拉 
E 吉尔吉斯 山 等 山脉 ,主要 的 河流 有 玛 纳 斯 河 、 伊 
Aam EW . 塔 拉 斯 河 等 。 本 文 以 河流 出 山口 为 划 
分 依据 提取 流域 边界 ,将 天 山北 坡 冰 川 区 划分 为 19 
个 子 流域 。 博 格 达 山北 坡 包括 开垦 河 .中 葛根 河 、 
白杨 河 等 流域 ,这 些 流域 范围 均 较 小 且 冰 川 分 布 相 
对 较 少 , 故 在 此 统一 合并 为 博 格 达 山 北 坡 流域 。 

天 山 山 脉 发 育 的 冰川 属于 亚 大 陆 型 冰川 或 大 
陆 型 冰川 ,形态 类 型 多 样 , 且 以 山谷 冰川 较 多 .5 
据 RGI(Randolph glacier inventory ) 6.0 KNE: 天 
EYE B UKJI] 7977 2& ,总 面积 约 
5059.67 km" , Kil Bil Bs v1 CU P E n DC, 
乌鲁木齐 河源 1 号 冰川 、Kara-Batkak 冰川 Ts. 
Tuyuksuyskiy 冰川 „Igli Tuyuksu 冰川 和 Golubin 冰川 
(图 1) 是 天 山北 坡 的 重点 监测 冰川 , 自 20 世 纪 50、 


人 


po Sul | " 


鲁 木 齐 河流 域 ;6: 头 屯 河 流域 ;7; 三 屯 河 流域 ; 


8: 呼 图 辟 河 流域 ;9: 塔 西河 流域 ;10: 玛 纳 斯 河流 域 ;11: 宁 家 河流 域 ;12: 人 金 沟 河流 域 ;13: ea 14: 伊 犁 河 -巴尔 喀什 湖 流域 
15: 艾 比 湖 流域 ;16: 阿 拉 湖 流域 ;17: 伊 塞 克 湖 流域 ;18 : 楚 河 流域 ;19: 塔 拉 斯 河流 域 。 
图 1 天 山北 坡 诸 河 流域 冰川 分 布 


Fig. 1 Distribution of glaciers in the sub-basins on the north slope of Tianshan Mountains 
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60 年 代 就 开始 有 物质 平衡 观测 记录 。Ts. Tuyuk- 
suyskiy 冰川 和 乌鲁木齐 河源 1 号 冰川 观测 序列 最 
长 ,是 国际 上 公认 的 全 球 重 点 监测 冰川 ,其 中 乌 鲁 
木 齐 河源 1 号 冰川 是 我 国 监测 时 间 最 长 ,观测 资料 
连续 性 最 好 的 一 条 冰川 。 


2 数据 与 方法 


2.1 数据 来 源 与 处 理 

2.1.1 Landsat OLI i£ & 3j 4$ — 本 文选 用 Landsat 
0LI 遥感 影像 (下 载 于 http://glovis.usgs.gov) 来 提取 
近期 天 山北 坡 冰 川 边界 ,数据 产品 为 Level 1T 级 , 空 
间 分 辩 率 为 30 m。 选 取 影 像 时 ,对 天 山北 坡 冰 川 区 
2016 年 消融 季 末 期 的 所 有 影像 进行 综合 对 比 , 选 取 
了 受 云 . 积 雪 和 山体 阴影 等 因素 干扰 最 小 的 影像 ， 
对 于 无 法 获取 消融 季 末 期 适用 影像 的 区 域 ,通过 选 
取 相 邻 年 份 的 影像 予以 补充 。 最 终 选 用 70 景 Land- 
sat OLL 影像 ( 表 1) ,影像 时 段 主要 集中 在 2014 一 
2017 ^F., HLEA 2015 年 为 主 , 故 将 本 文 资 料 时 期 统称 
为 2015 年 。 此 外 ,在 冰川 边界 解 译 之 前 ,通过 全 色 
波段 融合 将 影像 分 辩 率 提高 至 15 m, 以 更 好 地 满足 
冰川 边界 提取 的 精度 要 求 。 

2.1.2 DEM 数据 本 文选 用 美国 国家 航空 航天 局 
发 布 的 ASTER DEM V003(AST14DEM) 数 据 ( 获 取 
于 https;//search.earthdata.nasa.gov ) 用 于 分 析 天 山北 
坡 冰 川 表面 高 程 变 化 。 该 数据 是 使 用 可 见 光 和 近 
红外 传感器 获取 的 ASTER Level 1A 数据 的 立体 像 
对 生成 的 数据 产品 ,空间 分 辩 率 为 30 m, EE TREE 
范围 为 2000 一 2022 年 。 本 文 主要 选用 2000 年 和 
2020 年 消融 季 末期 . 受 云 覆盖 影响 较 小 的 数据 ,对 
于 无 法 获取 适用 数据 的 区 域 , 选 取 邻 近年 份 的 数据 
予以 补充 。 

2.1.3 冰川 编目 数据 本 文选 用 中 国 第 一 次 、 第 二 
次 冰川 编目 数据 ,以 及 RGI 6.0 冰 川 编目 数据 作为 天 
山北 坡 冰川 边界 提取 时 的 参考 数据 ,并 在 分 析 冰 川 
变化 时 使 用 了 该 数据 。 中 国 第 一 次 冰川 编目 以 20 
世纪 50 一 80 年 代 的 航 摄 地 形 图 和 航空 相片 为 主要 
数据 源 呈 ,其 中 天 山北 坡 航 片 的 主要 时 段 为 
1960s。 中 国 第 二 次 冰川 编目 (来 源 于 http://westdc. 
westgis.ac.cn ) 以 2006 一 2011 年 Landsat TM/ETM +3% 
感 影像 为 主要 数据 源 "…… ,天 山北 坡 所 用 影像 以 2007 
年 为 主 。RGI 6.0 冰 川 编 目 (来 源 于 http://www.glims. 
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org/RGI) 以 1999 年 及 之 后 的 Landsat TM/ETM+ 影 像 
为 主要 数据 源 "" ,中 国 部 分 均 为 第 二 次 冰川 编目 数 
据 ,该 数据 集中 天 山北 坡 国外 段 的 影像 时 段 以 1999 
年 为 主 。 在 分 析 各 子 流域 冰川 变化 时 ,以 其 所 用 影 
像 的 主要 时 段 来 代表 编目 的 资料 时 期 。 
2.1.4 Google Earth 影像 和 SAR 影像 本 文 使 用 
Google Earth 影像 和 SAR 影像 作为 冰川 边界 提取 的 
参考 影像 。Google Earth 可 查看 多 时 相 、 高 分 辨 率 的 
遥感 影像 , 旦 以 3D 形 式 显 示 了 高 原 、 山 地 等 地 形 ， 
而 能 够 清楚 识别 表 硬 覆盖 冰川 的 末端 、 侧 研 等 准确 
位 置 ,有 效 提高 解 译 精度 。SAR 因 具 有 穿 透 云 层 、 
全 天 候 、 部 分 穿 透 性 及 高 分 辨 紊 的 特征 ,对 难以 通 
过 光学 影像 识别 的 表 研 覆盖 冰川 ,参考 SAR 影像 干 
涉 失 相干 特性 进行 识别 。 
2.1.5 冰川 物质 平衡 观测 数据 本 文 在 分 析 乌 鲁 木 
齐 河 源 1 号 冰川 、Kara-Batkak 冰川 Ts. Tuyuksuyskiy 
冰川 Igli Tuyuksu 冰川 和 Golubin 冰川 长 期 物质 平 
衡 变 化 ,以 及 评估 采用 大 地 测量 学 方法 所 得 天 山北 
坡 冰 川 物质 平衡 估算 结果 的 准确 性 时 使 用 了 物质 
平衡 观测 数据 。 冰 川 物质 平衡 观测 数据 来 源 于 世 
界 冰 川 监测 服务 处 (Word Glacier Monitoring Ser- 
vice, WGMS) (https://wgms.ch/ggcb ) 和 相关 研究 "2 。 
2.2 研究 方法 
2.2.1 冰川 边界 提取 方法 利用 遥感 手段 提取 永川 
边界 的 方法 主要 分 为 目 视 解 译 和 计算 机 自动 分 类 2 
种 。 考 虑 到 天 山地 区 存在 大 量 的 表 研 履 盖 冰川 , 难 
以 通过 自动 分 类 方法 准确 识别 冰川 边界 ,因此 ,本 
文采 用 目 视 解 译 的 方法 提取 冰川 边界 。 对 Landsat 
0OLI 影 像 进行 假 彩色 合成 ,合成 后 的 影像 中 裸 冰 区 
与 非 冰 川 区 的 色彩 ,纹理 等 要 素 的 差异 十 分 显著 
(图 2), 因 此 ,对 于 裸 冰 区 而 言 ,可 直接 根据 冰川 区 
与 非 冰川 区 显著 的 色彩 差异 来 提取 冰川 边界 。 
而 对 于 冰川 的 表 厂 履 盖 区 ,结合 Google Earth 影 
像 根据 表 犬 覆盖 区 与 非 冰 川 区 的 色彩 和 纹理 差异 ， 
以 及 表 研 覆盖 区 的 地 形 与 水 文 特征 (如 冰川 末端 的 
拱 形 地 形 特征 、 冰 面 湖 的 分 布 等 ) 进 行 识别 。 对 于 
难以 通过 光学 影像 识别 的 表 研 冰川 ,采用 C 波 段 
Sentinel-1 Level 1 SLC 数据 和 AST14DEM 数据 进行 
差分 干涉 处 理 , 经 过 精 配 准 生 成 干涉 条 纹 ,滤波 消 
除 平 地 相位 URL DEM 模拟 相位 消除 地 形 相 位 的 过 
程 ,得 到 冰川 区 干涉 相干 系数 图 ,结合 光学 影像 与 
SAR 影像 干涉 失 相 干 特 性 来 进行 识别 (图 3)。 冰 川 
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表 1 Landsat 遥 感 影像 列表 


Tab.1 List ofLandsat remote sensing images 


影像 编号 ) 云 量 /% 影像 编号 po 云 量 /% 
LC81380302015243LGN00 2015-08-31 4.49 LC81470282016245LGN00 2016-09-01 0.04 
LC81380302016214LGN00 2016-08-01 6.21 LC81470292014239LGN00 2014-08-27 4.26 
LC81390302015202LGN00 2015-07-21 1.05 LC81470292016245LGN00 2016-09-01 0.14 
LC81390302019229LGN00 2019-08-17 6.22 LC81470292017263LGN00 2017-09-20 0.91 
LC81410302016267LGN00 2016-09-23 0.29 LC81470302014207LGN00 2014-07-26 1.88 
LC81410302017269LGN00 2017-09-26 0.03 LC81470302016245LGN00 2016-09-01 1.01 
LC81420302016210LGN00 2016-07-28 0.63 LC81470302017247LGN00 2017-09-04 1.64 
LC81420302016242LGN00 2016-08-29 4.20 LC81470312014207LGN00 2014-07-26 17.98 
LC81420302017244LGN00 2017-09-01 3.80 LC81470312015226LGN00 2015-08-14 13.77 
LC81420302017260LGN00 2017-09-17 0.48 LC81470312016245LGN00 2016-09-01 6.71 
LC81420302017276LGN00 2017-10-03 2.51 LC81480292015217LGN00 2015-08-05 2.00 
LC81430302014243LGN00 2014-08-31 0.99 LC81480292015233LGN00 2015-08-21 0.81 
LC81430302016217LGN00 2016-08-04 0.75 LC81480292016268LGN00 2016-09-24 2.26 
LC81430302017203LGN00 2017-07-22 1.26 LC81480312014214LGN00 2014-08-02 12.93 
LC81440302014202LGN00 2014-07-21 15.00 LC81480312015233LGN00 2015-08-21 7.09 
LC81440302014234LGN00 2014-08-22 1.46 LC81480312016252LGN00 2016-09-08 17.99 
LC81440302015285LGN00 2015-10-12 8.08 LC81480312017190LGN00 2017-07-09 5.78 
LC81440302017210LGN00 2017-07-29 19.59 LC81490302015208LGN00 2015-07-27 9.69 
LC81450292014241LGN00 2014-08-29 0.58 LC81490302015224LGN00 2015-08-12 1.79 
LC81450292015196LGN00 2015-07-15 1.34 LC81490302017261LGN00 2017-09-18 0.72 
LC81450292015212LGN00 2015-07-31 7.32 LC81490312015224LGN00 2015-08-12 0.99 
LC81450292016231LGN00 2016-08-18 9.59 LC81490312016259LGN00 2016-09-15 0.67 
LC81450292016263LGN00 2016-09-19 11.35 LC81490312017245LGN00 2017-09-02 1.54 
LC81450302014241 LGN00 2014-08-29 5.94 LC81500312015231LGN00 2015-08-19 0.50 
LC81450302015196LGN00 2015-07-15 0.50 LC81500312016266LGN00 2016-09-22 1.14 
LC81450302015212LGN00 2015-07-31 19.06 LC81510302014219LGN00 2014-08-07 0.10 
LC81450302016263LGN00 2016-09-19 3.43 LC81510302015254LGN00 2015-09-11 2.06 
LC81450312015196LGN00 2015-07-15 0.45 LC81510302017243LGN00 2017-08-31 0.12 
LC81460292015203LGN00 2015-07-22 6.10 LC81510312016241LGN00 2016-08-28 0.55 
LC81460292019246LGN00 2019-09-03 5.81 LC81520302017218LGN00 2017-08-06 0.18 
LC81460302014200LGN00 2014-07-19 1.03 LC81520312014242LGN00 2014-08-30 0.02 
LC81460302016222LGN00 2016-08-09 1.66 LC81530312014233LGN00 2014-08-21 0.07 
LC81460302018227LGN00 2018-08-15 6.74 LC81530312016239LGN00 2016-08-26 1.19 
LC81460312015203LGN00 2015-07-22 11.56 LC81530312018260LGN00 2018-09-17 3.15 
LC81460312016222LGN00 2016-08-09 2.50 LC81530312020234LGN00 2020-08-21 3.69 


边界 提取 后 ,对 结果 进行 了 进一步 的 检查 修正 , 计 2.2.2 冰川 面积 变化 计算 方法 考虑 到 不 同 冰 川 编 
算 了 相关 属性 信息 ,并 且 掏 除了 其 中 面积 <0.01 km ” 目 获 取 时 间 存 在 差异 ,为 确保 面积 变化 的 可 对 比 
的 对 象 。 性 ,本 文采 用 冰川 面积 变化 率 作为 衡量 面积 变化 的 


图 2 裸 冰 区 冰川 边界 识别 


Fig. 2 Identification of clean-ice boundaries 


Ey ” ,计算 公式 如 下 : 
MG 
x i 


式 中 ;AC 为 冰川 面积 变化 率 (%a'); AS; 为 第 i 时 
段 的 冰川 面积 变化 量 (km ); S; 为 第 ;时 段 内 起 始 年 
份 的 冰川 面积 (km ); AT, 为 第 ;时 段 的 时 间 间 隅 (a)。 


(a) 区 域 一 Landsat 遥 感 影像 


CS 冰川 边界 


(d) 区 域 二 Landsat 遥 感 影像 


CS 冰川 边界 


(b) 区 域 一 Google Earth 影像 
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2.2.3 物质 平衡 估算 方法 
(1) DEM 配 准 方法 
本 文 根 据 Nuth 等 提出 的 空间 匹配 误差 导致 的 
高 程 差 与 地 形 因 子 之 间 的 关系 来 进行 空间 位 置 偏 
差 的 校正 ,计算 公式 如 下 ， 
dh/tan (a)=aXcos(b-9)+e (2) 
c= dh/tan (a) (3) 
式 中 :qh 为 DEM 数 据 间 的 高 程 差 (m); e 和 a 分 别 
为 对 应 像 元 的 坡 向 和 坡度 (°);a 和 2 分 别 为 水 平 位 
移 的 偏 移 量 (m) 和 角度 (°); dh X DEM 数据 间 的 平 
均 高 程 差 (m)。 本 文 还 根据 高 程 差 的 5% 和 95% 分 
位 数 对 异常 值 进行 剔除 。 最 后 ,利用 非 冰 川 区 高 程 
差 残 差 与 最 大 曲率 之 间 的 线性 关系 ,对 冰川 区 的 高 
程 差 残 差 进行 校正 。 
(2) 物质 平衡 估算 方法 
利用 大 地 测量 学 方法 估算 冰川 物质 平衡 ,计算 
公式 如 下 : 


p n 
Bog X PARS, (4) 


8 


(c) 区 域 一 SAR 影 像 干涉 相干 系数 图 


图 3 表 太 区 冰川 边界 识别 


Fig.3 Identification of debris-covered glacier boundaries 
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式 中 :8 为 物质 平衡 (m w.e.-a); p 为 转换 密度 ; S, 
为 冰川 面积 (km) ;n 为 冰川 区 的 像 元 总 数 ; AR; 为 研 
究 时 段 内 像 元 点 的 高 程 差 (m); S, 为 单个 像 元 面积 
(km2)。 本 文选 用 Huss 研究 推荐 的 850+60 kg m^ ME 
为 转换 密度 。 
2.2.4 精度 评价 

(1) 冰川 面积 误差 评估 

本 文采 用 统计 冰川 轮廓 线 经 过 的 像 元 数量 的 
方法 来 评估 冰川 面积 的 误差 ,计算 公式 如 下 : 

E,-NxA (5) 

式 中 :及 为 冰川 面积 的 不 确定 性 (km) ; N 为 冰川 边 
界 ( 不 包括 用 于 冰川 分 割 的 山 状 线 或 分 冰 怜 ) 经 过 
的 像 元 数量 ;为 遥感 影像 像 元 大 小 (融合 后 的 Land- 
sat OLI 影 像 像 元 大 小 为 15 m); 

(2) 冰川 物质 平衡 误差 评估 

采用 非 冰 川 区 平均 高 程 差 和 平均 标准 误差 来 
评估 冰川 表面 高 程 变化 的 不 确定 性 ,计算 公式 
如 下 : 


SE=STDV/NK (6) 
o = NMED’ + SE? (7) 


式 中 :SE 为 平均 标准 误差 ;STDV 为 非 冰 川 区 DEM 
高 程 差 的 标准 差 ;K 为 空间 去 相关 处 理 (去 相关 距离 
设 定 为 600m) 后 的 像 元 个 数 ; o 为 DEM 数据 高 程 
变化 的 相对 误差 ;MED 为 非 冰 川 区 DEM 数据 的 平 
均 高 程 差 。 

人 研究 时 段 内 物质 平衡 的 不 确定 性 ( wy ) 计 算 公 
CP: 


Uy 一 


(Chappe (qu) 
式 中 : AR 为 冰川 区 的 平均 高 程 差 (m) ;为 研究 时 段 
长 度 (2000 一 2020 年 ); pv 和 p, 分别 为 水 的 密度 
(1000 kg*m”) 和 转换 密度 (850 kg m7) ; Ap 为 冰 密 
度 的 不 确定 性 (60 kgm”) o 


3 结果 与 分 析 


3.1 天 山北 坡 冰 川 面积 变化 
3.1.1 天 山北 坡 整 体 冰 川 面积 变化 2015 年 前 后 天 
山北 坡 诸 河流 域 共计 分 布 冰川 10061 条 ,总 面积 约 
4855.85+245.86 km。 其 中 ,天 山北 坡 国内 段 共计 分 
布 冰川 6097 条 ,总 面积 约 3045.09+146.49 km*。 中 
国 第 一 次 冰川 编 日 显示 ,1960s 天 山北 坡 国内 段 冰 
川 总 数量 为 5764 条 ,总 面积 约 4269.30 km ;中 国 第 
二 次 冰川 编目 显示 ,2007 年 该 区 域 冰 川 总 数量 为 
5223 条 ,总 面积 约 3299.28+181.68 km 。 通 过 对 比 
可 以 看 出 ,1960s 一 2015 年 ,天 山北 坡 国 内 段 冰 川 总 
面积 共 减 少 1224.21+146.49 km*( 约 28.67%+3.43%)， 
面积 萎缩 速率 为 0.52%…a'+0.06%a'( 表 2)。 其 中 ， 
1960s 一 2007 年 和 2007 一 2015 年 面积 蓉 缩 速率 分 别 
H 0.4896 - a^ 0.0996 - a^! FN 0.96% * a £0.8896 «a^ , K 
8j 32 psi o HE IE ER 5 D E S A ES A 

RGI 6.0 冰 川 编 目 显 示 ,1999 年 天 山北 坡 国 外 段 
共计 分 布 冰川 2754 条 ,总 面积 约 1760.39 km’, RGI 
6.0 在 天 山北 坡 存在 部 分 区 域 漏 画 的 问题 (图 4) ,本 
文 对 遗漏 的 冰川 予以 补充 ,经 统计 ,2015 年 天 山北 
坡 国外 段 共计 分 布 冰川 3964 条 ,总 面积 约 1810.77+ 
99.37 km 。 为 准确 分 析 该 区 域 冰 川 面积 变化 状况 ， 
将 本 文 提取 结果 中 天 山北 坡 国外 段 与 RGI 6.0 资 料 
能 够 完全 对 应 的 共计 2790 条 冰川 ( 约 占 天 山北 坡 国 
外 段 冰川 总 面积 的 76%) 进 行 对 比 ,结果 显示 ,这 些 
冰川 的 总 面积 从 1999 年 的 1505.03 km 萎缩 至 2015 
年 的 1369.08+75.78 km? , 共 减 少 了 135.95+75.78 km? 
( 约 9.03+5.04% ) , 面积 萎缩 速率 约 为 0.569%6 .at 
0.319% .am。 
3.1.2 天 山北 坡 各 子 流域 冰川 面积 变化 对 天 山北 
坡 各 子 流域 冰川 面积 变化 进行 统计 ,结果 显示 ， 
1960s—2015 年 天 山北 坡 国内 段 各 子 流域 冰川 面积 
均 呈 萎缩 趋势 ,其 中 头 地 河流 域 冰 川 面 积 萎 缩 速 率 


表 2 1960s 一 2015 年 天 山北 坡 国内 段 冰川 面积 变化 


Tab.2 Changes in glacier area on the north slope of Tianshan Mountains in China from 1960s to 2015 


时 段 变化 量 /km 相对 变化 率 /% 相对 变化 速率 /oa 
1960s 一 2007 年 -970.02+181.68 -22.72+4.26 -0.48+0.09 
2007—2015 年 -254.19+233.38 -7.70+7.07 -0.96+0.88 
1960s 一 2015 年 -1224.21+146.49 -28.67+3.43 -0.52+0.06 
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(a) 示例 一 Landsat 遥 感 影像 


图 例 
C3 RGI 6.0 冰 川 边界 
C3 2015 年 冰川 边界 


(b) 示例 二 Landsat 遥 感 影像 


CS RGI 6.0 冰 川 边界 
C3 2015 年 冰川 边界 
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图 4 2016 年 9 月 1 日 Landsat OLI 影 像 显示 的 RGI 6.0 资 料 与 本 文 提取 冰川 边界 对 比 
Fig. 4 Comparison between RGI 6.0 data and glacier boundaries extracted in this study 
by Landsat OLI image on September 1, 2016 


最 大 ,达到 1.17%…a'+0.11%…a ,其 次 为 乌鲁木齐 河 
流域 (1.16% - a^ ' x 0.1096 .aa 和 巴里 坤 湖 流 域 
(1.1696 * a! 0.0896 * a`) ,而 金 沟 河流 域 冰 川 面 积 
缩 速 率 最 小 (0.20%…a1!+0.07%a1)。 各 子 流域 冰川 
面积 普遍 呈 加 速 萎缩 趋势 , 头 地 河流 域 加 速 趋势 最 
明显 ,1972 一 2007 年 和 2007 一 2016 年 冰川 面积 萎缩 
率 分 别 为 0.99% - a^ 0.1796 - a^' FI 2.8396 .ai 
1.3396 -a'o JESP , 1999—2015 年 天 山北 坡 国外 段 
各 子 流域 冰川 面积 均 呈 萎缩 趋势 ,其 中 伊 塞 克 湖 流 
域 冰 川 面积 萎缩 速率 最 大 ,达到 0.87% .a'+t 
0.38% - a` ,而 塔 拉 斯 河流 域 冰 川 面积 萎缩 速率 相对 
较 小 (0.34%:a'+0.42%*a')。 

为 揭示 天 山北 坡 冰 川 面积 变化 的 空间 差异 ,将 
各 子 流域 近 十 几 年 来 的 冰川 面积 萎缩 速率 进行 对 
比 ,其 中 ,国内 和 国外 部 分 的 时 段 分 别 为 第 二 次 冰 
川 编目 和 RGI6.0 至 本 文 冰川 提取 时 段 , 在 此 统一 表 
述 为 2000s 一 2015 年 。 由 于 伊犁 河 -巴尔 喀什 湖 流 
域 婚 有 国内 部 分 也 有 国外 部 分 ,在 此 取 二 者 的 平均 
变化 速率 。 结 果 显 示 ( 图 5), 天 山北 坡 冰 川 面积 变 
化 具有 明显 的 空间 差异 , 东 段 各 子 流域 冰川 面积 
缩 速 率 普遍 较 西 段 更 大 ,其 中 头 屯 河流 域 获 缩 速率 
最 大 ,达到 约 2.83%a'+1.33%a', 其 次 为 乌鲁木齐 
河流 域 和 巴里 坤 湖 流域 ,面积 萎缩 速率 分 别 为 
2.5496 *a 1.2096 a^! 412.2796 a £1.1896*a^ 4. UKJI] 


面积 变化 的 空间 差异 与 冰川 规模 C RE ET 
素 有 关 ,气温 决定 冰川 消融 ,而 降水 量 决 定 冰 川 积 
累 。1960s 以 来 天 山北 坡 气温 显著 增加 , 东 段 各 子 
流域 升温 速率 相对 较 大 ,尤其 以 天 格 尔 山北 坡 气温 
升幅 最 为 显著 " ,该 区 域 冰 川 面积 蓉 缩 速率 也 最 大 ; 
降水 量 整体 呈 微 弱 的 增加 趋势 , 东 段 降水 量 明 显 少 
于 西 段 ,其 物质 补给 相应 地 较 少 ,冰川 面积 萎缩 速 
率 相对 较 大 。 此 外 ,冰川 面积 萎缩 速率 较 大 的 头 屯 
河和 乌鲁木齐 河流 域 冰 川 平 均 面 积分 别 仅 为 0.13+ 
0.01 km 和 0.16+0.01 km ,而 萎缩 速率 较 小 的 金 沟 河 
流域 冰川 平均 面积 约 0.82+0.03 km ,上 且 天 山北 坡 东 
段 冰川 规模 总 体 小 于 西 段 。 可 见 , 小 规模 冰川 对 升 
温 较 敏感 ,消融 更 强烈 ,可 能 是 造成 天 山北 坡 冰 川 
面积 变化 空间 差异 性 的 重要 原因 。 

3.2 天 山北 坡 冰 川 物质 平衡 变化 

3.2.1 天 山北 坡 监 测 冰 川 物质 平衡 变化 利用 冰川 
学 方法 观测 得 到 的 物质 平衡 结果 显示 (图 6), 过 去 几 十 
年 间 , 天 山北 坡 的 重点 监测 冰川 (图 1) 均 长 期 处 于 负 平 
衡 状态 。 乌 鲁 木 齐 河源 1 号 冰川 和 Ts. Tuyuksuyskiy 
冰川 长 期 连续 的 物质 平衡 观测 记录 显示 ,其 平均 物质 
平衡 分 别 为 -0.34 m w.e. +a 和 -0.38 m w.e."a'。 为 
验证 本 文采 用 大 地 测量 学 方法 所 得 天 山北 坡 冰 川 
物质 平衡 估算 结果 的 准确 性 ,将 本 文 估算 的 监测 冰 
川 物质 平衡 与 物质 平衡 观测 数据 及 前 人 研究 结 
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注 : 由 于 第 二 次 冰川 编目 数据 在 金 沟 河 流域 (12) 不 完整 ,所 以 其 面积 变化 未 比较 。 


图 5 2000s 一 2015 年 天 山北 坡 各 子 流域 冰川 面积 变化 


Fig.5 Changes in glacier area of sub-basins on the north slope of Tianshan Mountains during 2000s—2015 
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图 6 近 几 十 年 来 天 山北 坡 重点 监测 
冰川 累积 物质 平衡 变化 


Fig.6 Changes in cumulative mass balance of monitoring 


glaciers on the north slope of Tianshan 


Mountains in recent decades 


进行 对 比 ( 表 3)。 结 果 显 示 , 不 同方 法 获取 的 冰川 
物质 平衡 具有 较 好 的 一 致 性 ,表明 本 文物 质 平衡 估 
算 结果 具有 一 定 的 可 信和 度 ,而 不 同 结果 之 间 存 在 一 
定 差异 的 原因 可 能 与 研究 时 段 不 同 . DEM 数据 精度 
有 限 等 因素 有 关 。 

3.2.2 天 山北 坡 整 体 及 各 子 流域 冰川 物质 平衡 变 
化 基于 2000 年 和 2020 年 左右 的 AST14DEM 数据 ， 
采用 大 地 测量 学 方法 对 整个 天 山北 坡 冰 川 物 质 平 
衡 变化 进行 估算 ,结果 显示 ,2000 一 2020 年 该 区 域 
冰川 表面 高 程 明显 下 降 , 其 变化 速率 为 -0.57+ 
0.01 ma ,冰川 物 质 呈 持续 亏损 态 ,物质 亏损 速率 


达到 0.39+0.04 m w.e.…a'。 对 天 山北 坡 各 子 流域 冰 
川 物 质 平衡 变化 进行 统计 ,结果 表明 ,各 子 流域 冰 
川 物质 均 呈 亏损 状态 ,但 不 同 流域 的 物质 平衡 变化 
存在 一 定 的 空间 差异 性 (图 7)。 其 中 ,位 于 哈 尔 里 
克 山 北 坡 的 伊 吾 河流 域 冰川 物质 亏损 速率 最 大 , 达 
到 0.68+0.06 m w.e…a ,其 次 为 头 屯 河 流域 (0.67=+ 
0.07 m w.e. * a) ,而 宁 家 河流 域 物质 亏损 速率 最 
小 。 按 照 RGI 资 料 中 天 山东 、 西 段 的 划分 依据 统计 
显示 ,天 山北 坡 东 、 西 段 物 质 平 衡 估算 结果 分 别 
为 -0.42+0.04 m w.e. a 和 -0.37+0.04 m w.e."a', 可 
见 东 段 物 质 亏损 速率 相对 较 大 。 

此 外 ,本文 还 将 冰川 物质 平衡 估算 结果 与 前 
人 研究 结果 进行 了 对 比 。Brun 等 基于 ASTER 
立体 像 对 估算 得 到 2000—2016 年 整个 天 山地 区 
冰川 物质 平衡 为 -0.28+0.20 m w.e.- a! ,其 中 东 Lu 
段 物质 平衡 分 别 为 -0.40+0.20 m w.e. ar 和 -0.20+ 
0.08 mw.e.…a1。 该 研究 与 本 文 在 山系 分 段 上 存在 一 
定 差异 ,但 其 划分 的 天 山东 段 大 部 分 位 于 本 文 天 山 
北 坡 范围 内 ,二 者 在 东 段 的 物质 平衡 估算 结果 也 较 
为 接近 。 陈 安安 2 利用 SRTM DEM fI ASTER 立体 
像 对 数据 研究 发 现 ,1999 一 2015s 位 于 天 山北 坡 的 
伊犁 河流 域 冰 川 物 质 平衡 为 -0.42+0.16 m w.e.*a^', 
与 本 文 估算 结果 (-0.39+0.04 m w.e.…a') 基 本 一 致 。 
此 外 ,关于 天 山地 区 物质 平衡 空间 差异 性 的 诸多 研 
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表 3 天 山北 坡 重 点 监测 冰川 物质 平衡 结果 对 比 
Tab. 3 Comparison of glacier mass balance for monitoring glaciers on the north slope of Tianshan Mountains 
监测 冰川 时 段 方法 物质 平衡 /m w.e. +a” 数据 来 源 
乌鲁木齐 河源 1 号 冰川 2002 一 2017 年 冰川 学 方法 -0.66 WGMS 
2012 一 2018 年 大 地 测量 学 方法 -1.13+0.18" [23] 
2002 一 2017 年 大 地 测量 学 方法 一 0.48+0.04 本 研究 
Ts. Tuyuksuyskiy 冰川 2003—2018 4E 冰川 学 方法 -0.50 WGMS 
1998 一 2016 年 大 地 测量 学 方法 -0.35+0.18 [24] 
2003—2018 年 大 地 测量 学 方法 -0.51+0.04 本 研究 
Igli Tuyuksu 冰 川 1998 一 2016 年 大 地 测量 学 方法 -0.37+0.16 [24] 
2003 一 2018 年 大 地 测量 学 方法 -0.37+0.04 本 研究 
Golubin 冰 川 2010 一 2017 年 冰川 学 方法 -0.38 WGMS 
2000 一 2019 年 模型 估算 方法 -0.30 [25] 
2003 一 2017 年 大 地 测量 学 方法 —0.38+0.03 本 研究 
Kara-Batkak 冰川 2013 一 2020 年 冰川 学 方法 -0.70 WGMS 
2000 一 2018 年 模型 估算 方法 一 0.54+0.08 [26] 
2004—2020 年 大 地 测量 学 方法 -0.45+0.05 本 研究 
注 :* 表 示 冰 舌 区 物质 平衡 变化 ; WGMS 为 世界 冰川 监测 服务 处 。 
河流 “一 流域 边界 冰川 物质 平衡 /m w.e.a7 
富 湖泊 ”一 天 山北 坡 界 线 | 0- 100km 
es 冰川 7034 
图 7 2000 一 2020 年 天 山北 坡 各 子 流域 冰川 物质 平衡 变化 
Fig.7 Glacier mass balance of sub-basins on the north slope of Tianshan Mountains during 2000—2020 
究 结 果 均 显 示 , 近 十 几 年 间 天 山东 段 冰川 物质 亏损 。 冰川 物质 平衡 变化 特征 ,得 到 以 下 主要 结论 : 
速率 略 大 于 西 段 ”", 与 本 文 研究 结果 基本 一 致 。 (1) 2015 年 前 后 天 山北 坡 诸 河流 域 共计 分 布 冰 


论 


川 10061 条 ,总 面积 约 4855.85+245.86 km", 1960s— 


4 结 2015 年 天 山北 坡 国内 段 冰川 总 面积 共 减 少 1224.21+ 
146.49 km ( 约 28.67%+3.43%) ,面积 萎缩 速率 为 


本 文 基于 Landsat OLI 遥 感 影 像 .Google Earth 影 
像 和 SAR 影 像 ,提取 天 山北 坡 诸 河流 域 近 期 冰川 边 
界 ,结合 中 国 第 一 次 .第 二 次 冰川 编目 和 RGI 6.0 vk 
川 编目 数据 ,分 析 天 山北 坡 近 60 a 来 冰川 面积 变化 
特征 ,并 基于 AST14DEM 数据 分 析 近 20 a 来 该 区 域 


0.5296 - a^ 0.0696 - a^, ELE 4E 2c 5E JL E 2S Apo A 
(0.9696 - a 0.8896 * a! ) ; 1999—2015 年 天 山北 坡 
国外 段 冰川 面积 萎缩 速率 约 为 0.36%.a :+0.31%.a :。 
天 山北 坡 冰川 面积 变化 具有 明显 的 空间 差异 ， 
2000s 一 2015 年 东 有 段 冰 川 面积 萎缩 速率 相对 较 大 ， 
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速率 最 大 (2.83%a!'+1.33%a!)， 


而 西 段 各 子 流域 冰川 面积 萎缩 速率 相对 较 小 。 


(2) 2000 一 2020 年 天 山北 坡 冰 川 表面 高 程 明显 


下 降 , 其 变化 速率 为 -0.57+0.01 m- a^ ,该 区 域 冰 
川 物质 呈 持 续 亏 损 态 ,物质 亏损 速率 达到 0.39+ 


0.04 m w.e."a so 
的 


天 山北 坡 冰 川 物质 亏损 存在 一 定 


空间 差异 性 , 东 段 各 子 流 域 物质 亏损 速率 相对 较 


大 (0.42+0.04 m w.e. -a` ) ,其 中 伊 吾 河流 域 物质 亏损 
速率 最 大 (0.68+0.06 m w.e.*a^) ,而 西 段 各 子 流域 冰 
川 物质 亏损 速率 相对 较 小 。 
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Abstract: Understanding glacier changes on the north slope of the Tianshan Mountains on the basin scale is criti- 
cal for the rational development and use of water resources in downstream oases. The recent glacier boundaries 
on the north slope of the Tianshan Mountains were obtained based on multisource remote sensing images. Com- 
bined with the glacier inventory released earlier and ASTER DEM, changes in the glacier area and mass balance 
in this region were analyzed. The results revealed that: (1) A total of 10061 glaciers were distributed on the north 
slope of the Tianshan Mountains around 2015, with a total area of 4855.85::245.86 km'. The glacier area on the 
north slope of the Tianshan Mountains in China decreased at a rate of 0.5296-a '+0.06% .a ' from 1960s to 2015. 
Furthermore, an accelerated decreasing trend was observed (0.9696 .a '+0.88% -a '). The glacier area in the for- 
eign part of the north slope of Tianshan Mountains decreased at a rate of 0.5696 .a '+0.31% .a ' from 1999 to 
2015. (2) The change rate of glacier surface elevation on the north slope of Tianshan Mountains from 2000 to 
2020 was -0.57+0.01 m-a ', and the glacier mass balance was -0.39+0.04 m w.e.-a '. (3) An obvious spatial dif- 
ference exists in the changes of glacier area and mass balance in the eastern and western parts of the north slope 
of Tianshan Mountains, and high decreasing rates of glacier area and mass balance were observed in the eastern part. 


Key words: north slope of Tianshan Mountains; glacier change; mass balance; remote sensing 


